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Abstract 

We provide an elementary proof of the quantum adiabatic theorem. 

1 Introduction 

The model of adiabatic quantum computation is a new paradigm for designing quantum algorithms, pro- 
posed by Farhi et al. [5|. It was recently established that this model is polynomially equivalent to the 
standard model of quantum circuits fT2l [Tll. Nevertheless, this model provides a completely different way of 
constructing quantum algorithms and reasoning about them. Therefore, it is seen as a promising approach 
for the discovery of substantially new quantum algorithms. 

Farhi et al. [ 4 1 have numerically studied adiabatic quantum algorithms on random instances of problems 
such as SAT and finding fc-cliques, with promising results for small instances. Rigorous results are, however, 
quite scarce. Van Dam et al. [ 12 1 and Reichardt [ 10] have constructed examples of SAT formulas for which a 
certain natural adiabatic algorithm performs poorly. This suggests that, at least for certain cases, the quantum 
adiabatic algorithms are not much stronger than the classical method of simulated annealing. On the other 
hand, there are instances on which classical simulated annealing fails but the adiabatic quantum algorithm 
succeeds [3]. A rigorous analysis of adiabatic algorithms in the general case appears to be difficult. 

The model of adiabatic quantum computation is based on a theorem known as the quantum adiabatic 
theorem [8]. Informally, this theorem says that if we take a quantum system whose Hamiltonian slowly 
changes from H\ to H2, then, under certain conditions on Hi and H2, the ground (lowest energy) state 
of Hi gets transformed to the ground state of H2. This is used to construct adiabatic algorithms for op- 
timization problems, in the following way. We take a Hamiltonian Hi whose ground state \%Jji) we know 
and a Hamiltonian H2 whose ground state corresponds to the solution of our optimization problem. Then, 
starting a quantum system in the state and slowly changing the Hamiltonian from Hi to H2 will solve 
our optimization problem. 

The adiabatic theorem has several proofs in the physics literature (see, e.g., |6l|2]|8]])- However, these 
proofs are rather involved and seem to give very little intuition. 
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Even the correctness of the adiabatic theorem has recently been questioned. Marzlin and Sanders Q 
and Tong et al. [11] have given counterexamples to some variants of the adiabatic theorem that were widely 
assumed to be true. Reichardt [ 10 1 recently claimed that none of the proofs examined by him contains a 
rigorous analysis of the convergence time. He includes in his paper another proof of the adiabatic theorem 
where he addresses this issue. His proof, however, follows the structure of previous proofs and does not 
seem to be more intuitive. 

In this paper, we give a new proof of the adiabatic theorem. Unlike all previous proofs, our proof is 
elementary and should be much more accessible to computer scientists. Moreover, we believe that our proof 
gives a good insight into why the adiabatic theorem holds: essentially, the proof shows that the error in 
the adiabatic evolution can be written as a certain geometric sum and that if the evolution is performed 
slow enough, this geometric sum almost completely cancels out. We hope that our proof will lead to new 
adiabatic algorithms, a better understanding of existing algorithms and contribute to settling the controversy 
about the correctness of the adiabatic theorem. 

2 Overview 
2.1 Main result 

Let H(s), < s < 1, be a Hamiltonian dependent on a parameter s. We refer to H as a time dependent 
Hamiltonian. We think of H(0) as the initial Hamiltonian and of H(l) as the final Hamiltonian. For a time 
dependent Hamiltonian H, we use the notation \\H\\ to denote max sg [ ,i] l!^( s )ll where || • || is the usual 
operator norm. We use a similar notation to denote the maximum norm (or absolute value) of other time 
dependent expressions. 

Let ty(s) be an eigenstate of H(s) with eigenvalue j(s) (in most applications, ^(s) is chosen to be the 
ground state of H(s)). When we say that we apply the adiabatic evolution given by H and \& for time T 
we mean that we initialize a system in the state ^(O) and then apply the continuously varying Hamiltonian 
H(t/T) for times t € [0, T]. We expect the final state of the system to be close to ^(1). Our main result is 

Theorem 2.1 Let H{s), < s < 1, be a time dependent Hamiltonian, let ^>{s) be one of its eigenstates, 
and let 7(s) be the corresponding eigenvalue. Assume that for any s € [0, 1], all other eigenvalues of H(s) 
are either smaller than j(s) — A or larger than 7(5) + A (i.e., there is a spectral gap of X around j(s)). 
Consider the adiabatic evolution given by H and ^ applied for time T. Then, the following condition is 
enough to guarantee that the final state is at distance at most 5 from ^(l)." 



In particular, this implies that as long as H has a 1/poly spectral gap around 7, we can reach a state that 
is at most 1/poly away from ^(l) in polynomial time. We remark that it might be possible to improve the 
dependence on A to A 3 or even A 2 . 

2.2 Overview of our proof 

The main part of our paper is concerned with the special case of Theorem 12. II in which 7(5) = for all 
< s < 1 (see Figure 0. The proof of this special case contains most of the important ideas and allows us 
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to avoid a few technical issues. Later, in Section|4j we complete the proof of Theorem l2.1l bv showing how 
the general case reduces to this special case. In this overview, we concentrate on the proof of the special 
case. In order to emphasize the high level structure of the proof, some details are omitted. 




Figure 1 : Special case of constant eigenvalue with a gap of A on both sides 

We start by discretizing the adiabatic evolution. Namely, we replace H(s) by a sequence of (time- 
independent) Hamiltonians Hq, . . . ,Hl-i each of which is applied for a small interval of time e = |J. 
Equivalently, we are applying the sequence of unitary transformations Uq = e lsHo , U\ = e ieHl Ul-i = 
e ieH L-i_ L e t g. — \^(j/i^ be the corresponding discretization of \P(s). Our goal has now become the 
following: show that the unitary transformation Ul-\ ■ ■ ■ Uq transforms go into a state close to g^- 




Figure 2: gj and Wj 

To show that, we consider a sequence w\ , . . . , wl where Wj + ± is defined as the projection of gj — gj + i 
to the subspace orthogonal to (see Figure |2j. We will show that 

u j9j = 9j = 9j+i + Wj+l + 0(1/ L 2 ) 
where 0(1/ L 2 ) denotes a vector of norm 0(1/ L 2 ). Then, the final state is 

L 

U L -i ■■■U g = g L + J2 Ul-i • • • U jWj + 0(1/ L). 



Showing that this expression is close to gi is equivalent to showing that the norm of Y^j=i Ul-i ' • • UjWj 
is small. We show this by proving that all but at most a small fraction of this sum cancels out. To show 
cancellations, we split the sum into smaller groups of A terms each. We then show that the norm of each 
group is at most <5A/(2L). 

For simplicity, consider the first group 

A 

Y, r > i •••''./''> (i) 

Since all terms start with the unitary Ul-i • • • Ua, the norm of (Q is the same as the norm of 

A 
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Next, we show that the wj's and the Ufs change relatively slowly. More precisely, we show that if A is 
sufficiently small compared to L, then we can make the following approximations: 

• replace all Wj by w\ ; 

• replace all Uj by U\. 

Thus, we obtain that the norm of Q is closely approximated by the norm of 

A-l 

Y, u i w i- ( 2 ) 

We now arrive at the heart of the proof. Express w\ as a sum of eigenvectors of U\,w\ = Ylt=i a k ( t ) k- 
Let Afc be the eigenvalue of H\ corresponding to the eigenvector cpk- Then, the above sum can be written as 

d A-l 
k=l j=0 

Recall that w\ is orthogonal to g\, and that g\ has eigenvalue in H\. Since we assumed H\ has a spectral 
gap of A, all the A^'s in the above sum (ignoring terms for which = 0) are at least A in absolute value. 
Hence, if we pick A large enough compared to j-, then most of the sum 

T,f=o eijXk£ cancels out, giving 

the desired result. These cancellations in the geometric sum are the essential reason why adiabatic evolution 
works. 

In the next sections, we make these arguments precise. 

3 Proof of a Special Case 

In this section we prove Theorem l2.1l for the special case in which the eigenvalue of the eigenvector that we 
follow is always (see Figure This case already captures the essential ideas in our proof. In Section |3] 
we will show how to reduce the general case to this special case. 

Before we begin, we need to address a minor technical issue. Given some Hamiltonian with an eigen- 
vector ^(s), we would like to say that the adiabatic evolution closely follows ^(s) in the h norm. However, 
notice that the phase of ^(s) is arbitrary. So, for example, &(s) = e ts ^(s) is an equally good eigenvector 
and clearly, the adiabatic evolution cannot be close to both ^(s) and to <£(s) in the I2 norm as the distance 
between them is large. A possible solution is to use a distance measure that is insensitive to global phase. 
We choose to take a different approach: we find a way to set the phase of ^(s) so that the adiabatic evolution 
closely follows ^(s) in the I2 norm. As it turns out, the correct way to choose the phase is such that for all 
s, (^'(s), ^(s)} = 0. In the next claim, we show that this is possible. It can be seen that for any unit vector 
^(s), this inner product is a complex number that has zero real part. Intuitively, it indicates the speed by 
which the phase of the vector rotates. 

Claim 3.1 Let ^(s), < s < 1, be a time-dependent unit vector in some Hilbert space, such that ^(s) is a 
differentiable function of s. Then, there is another time-dependent unit vector &(s) that is identical to ^(s) 
up to phase such that ($'(s), =0 and $(0) = ^(O). 
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Proof: Write $(s) = e il3 ^^(s) for some function j3 : [0,1] — ► R. Taking derivative, we have 

$'( s ) = e^)^'(s) + ie i/3(s) j 9'(s)*(s). 

Taking the inner product with Q(s), we obtain 

($'( S ),1>( S )} = (^( S ),^( S )}+i/3 / ( S ) 

since (^(s), ^(s)} = 1- In order to make this expression zero, we choose 

0(s) = f i(*'(t),*(t))dt. 







We also need the following technical lemma. Essentially, it says that if H changes slowly, then so does 
ty. We postpone its proof to Subsection l3.ll 

Lemma 3.2 Let H(s) be a time dependent Hamiltonian and let ^f(s) be an eigenvector with eigenvalue 0. 
Assume that the phase of the eigenvector is chosen such that (^'(s), ^(s)) = for all < s < 1. Moreover, 
assume that all other eigenvalues are at least A in absolute value. Then, 



and 

„ ..„ \\H"\\ 3\\H'\\ 2 
Kp" < I] !i + -J 1]_ 

11 A A 2 

The following is the main result of this section. 

Lemma 3.3 Let H(s), < s < 1, be a time dependent Hamiltonian and let 5 > be any constant. Let ^f(s) 
be an eigenvector whose eigenvalue is such that (^f'(s), ^f(s)) =0 and assume that all other eigenvalues 
are at least A in absolute value. Let T denote the time along which we apply the Hamiltonian. Then, the 
following condition is enough to guarantee that an adiabatic evolution starting from ^(O) is within distance 
6ofV(l): 

1000 f||#l 3 ll-ff'H • ||#"|| 1 
T ^' maX {^' A^ J 

Proof: For ease of presentation, we discretize time into infinitesimally small units of size 1/L. One should 
think of L as a quantity going to infinity while all other quantities remain constant. We use the 0() notation 
to describe the asymptotic behavior of an expression as a function of L; all other quantities are regarded as 
constants, e.g., \\H'f/L = 0(1/ L). 

We start by discretizing the adiabatic evolution. Let 

Uj := e i-T/L-H(j/L) 

be the unitary obtained by applying H(j/L) for T/L time units. Then the adiabatic evolution is closely 
approximated by the unitary 

U L -i ■■■U q . 
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The error in this approximation goes down to with L and can therefore be ignored (see, e.g., [ 12 1). 

Next, let gj = ty(j/L) be the discretized eigenvectors. In other words, gj is the eigenstate with eigen- 
value of H(j/L). The adiabatic evolution starts with go = \&(0). Notice that Ujgj = gj. Define 

W 3+l = P g± (9j-9j+l) 

where P g ± is the projection on the subspace orthogonal to gj + % (see Figure©. By taking the Taylor series 
of \F about (j + 1)/L, we obtain 

By applying P ± to both sides of the equality, we get 



9j+l 



Wj+i = _ n<j + m) +om% (3) 

where the ^'((j + 1)/L) term remains unchanged because we chose (^'(s), \&(s)) = for all s and, in 
particular, for s = (j + 1)/L. By combining the two equations above, we obtain 

9j = 9j+i + w j+1 + 0{1/L 2 ). 

Therefore, we can write 

U L -i ■ ■ ■ UiU g = f7 L _i • • • Uxgo 



2\ 



U L -. 1 ---U 1 (g 1 +w 1 ) + 0(l/L 



L 

gL + Y, u ^--- U i w 3 + °( 1 / L )- 

3=1 



Our goal is to show that the above is very close to g^ = ^(1). The term 0(1/ L) is negligible since L goes 
to infinity. Therefore, it is enough to show in the following that 



Yl Ul ~ 1 ' ' ' U 3 W 3 
3=1 



< 6. (4) 



First, according to Lemma U!2l 

HII < (5) 

Hence, if we try to bound the left side of © by a straightforward application of the triangle inequality we 
obtain a bound of SL -r^. In the remainder of the proof we will show how to improve this bound to 5. 

Define A as \(8/5)L\\H'\\/(TX 2 )]. Notice that A = 0(L). We start by partitioning the sum in © to 
sections containing A terms each. Namely, © will follow by showing that for any k, 

fc+A-l 6A 



j=k 



< 

~ L 
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For simplicity, let us consider the case k = 1; essentially the same proof works for any k. So, in the 
following we will show that 

II A 5A 

3=1 

Since Ul-i ■ ■ ■ Ua are unitary and are applied to every component of this sum, this is equivalent to 



3=1 



< 



SA 
L ' 



(6) 



< 



SA 
2L' 



(7) 



Later, we will show that 

A A-1 

|| U A -i ■ ■ • U jWj - Y U(W! 

j=l j=0 

Assuming ©, we can now complete the proof of the theorem. Let g be any eigenvector of H(l/L) such that 
g ^ g\ and let a denote its eigenvalue. The corresponding eigenvalue of g in U\ is e iaT / L . Since g ^ g±, 
\a\ > X and we can write 



A-1 A-1 

Yu{g =\Y^ T/L 



j=0 



j=0 



|iaAT/L_ 1 | 4L 4L 
< < 

~ T a ~ T\ 



(8) 



where we used that \e ie - 1| > |0|/2for any small enough 6. These cancellations in the geometric sum are 
at the heart of the adiabatic theorem. Using ©, we obtain that 



A-1 

E ^ 

3=0 



4L \\H'\\ 5A 
~ TX ' XL ~ 2L 



where the first inequality follows by writing w\ in the basis of eigenvectors of H(l/L) and recalling that 
is orthogonal to g\. Combined with Q, this proves © and completes the proof of the theorem. 
It remains to prove (Q. We will prove it in two steps. First, we will show that we can replace all Wj's 

with w\ and later we will show that we can replace all Uj's with U\. 

Lemma 3.4 For all j, k, 



\w j+k -Wj\\ < 



L 2 



\H"\\ \\H 



A 2 



0(1/L 2 



Proof: Using Eq. @, 



*>j+k ~ w i = i' + k )l L ) ~ *'U/L)) + 0(l/L 2 ). 

By the mean value theorem and the second claim in Lemma l3~2l the norm of the above is at most 



A||*1 + 0(1/L2)<A 



l-fT'll 11^ 



ni2 



A 2 



+ 0(1/L 2 ). 
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By the above lemma and the triangle inequality, we obtain 



A A A 

Ua-i ■ ■ ■ UjWj - ^ Ua-i ■ ■ ■ UjWi < ^ II^a-i • • • UjWj - Ua-i • • • UjW± \ 
j=x j=i j=i 

A 

= SH-Wlll 

i=i 



A 2 / \\H'- 



/||2 



* 1? ^T + 3i A 2 



+ 0(1/L) 



< 



5A 
IE 



where the last inequality is by our choice of T. 

In order to complete the proof, it is enough to show (notice that Ylf=o ^1 = S_/=i ^1 J ) : 



A A 

£ C/a-i • • • C/^i - ^ C/f" V 



< 



SA 
4L' 



Lemma 3.5 For all j, 



\(' r n-i: i \\<^p-+0(l/L 3 ) 



L 2 



Proof: Let J = H((j + 1)/L) — H(j/L). Then, using the Trotter formula [9], we can write 

Uj+1 = e i-T/L-H(( j+ l)/L) = e i.T/L.HU/L) e i.T/L.J + (1/L 3 ) 

where we used ||J|| = 0(1/ L). Then, 

\\U j+1 - Uj\\ = \\e l - T ' L - J - I\\ + 0(l/L 3 ) = j\\J\\ + 0(1/L 3 ) < + 0(1/L 3 ). 



(9) 



By combining this lemma with the triangle inequality, we obtain that for all k, 

kTW H'W 

||t4-£/i||<^-^+0(A;/L 3 ). 

We now prove Equation [9] using a sequence of A — 1 triangle inequalities, as illustrated in the following 
diagram: 

Ua-i ■■■U 3 U 2 U 1 Ua-i ■■■U 3 U 1 U 1 Ua-i ■■■U 1 U 1 U 1 U 1 ---U 1 U l U x 

Ua-i ■■■U 3 U 2 U A -i ■■■U 3 U 1 U A -i ■■■U 1 U 1 U 1 ---U 1 U 1 

Ua-i ■■■U 3 -h- Ua-i ■ • • U 3 -» Ua-i ■■■Ui -> . . . U x ■ ■ ■ U x 



Ua-i Ua-i U A -i C/i 
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That is, we use the triangle inequality to bound the left side of © by the sum of A — 1 terms where the 
k'th term is given by (notice that all terms not containing cancel and that the unitaries Ua-i, • • • , f^+i 
appear in all remaining terms and can therefore be ignored): 



3=1 3=1 



(Uk-Ui^Ut'wi 
3=1 

k 

< iif/ fc -t/iii • lE^r 1 ^ 



< 



L 2 
4fc||ff'|| 2 
L 2 A 2 



TX XL 



+ 0(k/L 3 



where the last inequality follows from (J3J and an argument similar to the one used after ©. Summing over 
k = 1, . . . , A — 1 we obtain that the left side of ^ can be upper bounded by: 

wa 2 \\h'\\ 2 



X 2 L 2 



< 



5A 
4L 



by our choice of T. ■ 

3.1 Proof of Lemma EI] 

The equality 

H{s)^{s) = 
holds for all < s < 1. Take the derivative according to s, 

H'(s)V(s)+H(s)y'(s) = 0. (10) 

and by taking the norm we obtain, 

\\H(s)y'(s)\\ = \\H'(s)^(s)\\ < \\H'\\. 

On the other hand, 

\\H( S )V(s)\\>X\\P ns) ^'( S )\\ = X\\V( S )\\ 

where we used the fact that all other eigenvalues of H(s) are at least A in absolute value and P^t 8 \± denotes 
the projection on the space orthogonal to ^(s). By combining the two inequalities, we obtain the first claim. 
For the second claim, let us consider the derivative of Equation fTOl 

#"(s)¥(s) + 2 • H'(s) • + H(s)V"(s) = 

and by taking the norm we obtain 



\H(s)V"(s)\\ < \\H"(s)*(s)\\ + 2 • \\H'(s) • *'(s)|| < \\H"\\ + 2 



\H'\ 



X 



On the other hand we have, 



\H{sW(s)\\>X\\P ns) ^"{s 
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from which we obtain 

\\P* {s) ±*"(s)\\ < M + (11) 
Now, by taking the derivative of (^(s), \&'(s)) = 0, 

and hence, 

K*( S ),*"( S ))| = ||*'( S )|| 2 . 

We complete the proof by combining the last equality with ( fTTT) : 

< \ (^( s ), + ||-FW( s )-l*"0)|| (using the triangle inequality) 



<||tt'( s )|f + JL_i> +2 a , 
\\H'\\ 2 

< — h 3 — -5 — (using the first claim). 

A A^ 

4 Reducing to a special case 

Lemma 4.1 Let H(s) be a time dependent Hamiltonian and let ^(s) be an eigenvector with eigenvalue 
j(s). Then, for any < s < 1, 

yw < ii^'ii 

and 

7"(s) < \\H"\\+4\\H'\\ 2 /X 

Before proving this lemma, let us see why it implies the main theorem. 

Proof of Theorem l2.H Define the Hamiltonian H(s) = H(s) —7(a)/. Since H and H differ by a multiple 
of the identity, they both describe the same adiabatic evolution up to some global phase. Moreover, by 
Lemma l4.ll 

\\H'\\ < \\H'\\ + IVI < 2||i7'|| 

and 

||#"|| < \\H"\\ + | 7 "| < 2\\H"\\+A\\H'\\ 2 /X. 
Hence, according to Lemma l331 it is enough to choose T to be at least 

1000 fll^'ll 3 \\H'\\ ■ \\H"\\ 1 10 5 (WH'P \\H'\\-\\H"\ 

max < — — ; — , > < —TT- ■ max 



5 2 A 4 ' A 3 - <5 2 A 4 ' A 3 



Proof of Lemma ETTl Using the Taylor expansion, we can write: 

j(ds) = 7(0) + ds7'(0) + ^sV(O) + 0(ds 3 ) = ds 7 '(0) + -ds 2 j"(0) + 0(ds 3 ) 
H{ds) = H(0)+dsH'{0) + ^ds 2 H"(0) + O(ds 3 ) 
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We know that 

H(ds)\^(ds)) =j(ds)\y(ds)) 
and by multiplying with (^(0)| we obtain: 

(*(0) | JET(ds) | * (d») ) = 7(<fa)(*(0)[*(ds)) 

With the above equalities, this simplifies to 

(V(0)\dsH'(p) + l/2ds 2 H"(0) + 0(ds 3 )\* (ds)) = (ds 7 '(0) + l/2ds 2 7 "(0) + O(ds 3 )){^(0)\^{ds)) 

and we can divide by ds: 

(V(0)\H\0) + l/2dsH"(0) + 0(ds 2 )\^(ds)) = ( 7 '(0) + l/2dsj"(0) + 0(ds 2 )){^> (0)|*(ds)) 

Since (^'(0)| v I / ((fs)) = 1 — 0(<is 2 ) we can hide all error terms inside the 0(ds 2 ): 

{$(0)\H\0) + l/2dsF"(0)|^(ds)} = 7 '(0) + l/2ds 7 "(0) + 0(ds 2 ) 

Now, ^(ds) = f(0) + dstf'(O) + Of/is 2 ): 

(*(0)|H'(0) + l/2dsH"(0)|*(0)) + (*(0)\dsH'(Q)\y'(Q)) = 7 '(0) + 1/2^/(0) + 0{ds 2 ) 

So by equating the coefficients of the polynomials, we obtain 

7 '(0) = <*(0)|IT'(0)|*(0)) < Hil'll 

and 

7 "(0) = (tf(0)(iZ"(0)|¥(0)) + 2<^(0)|#'(0)|^'(0)} < \\H"\\+2\\H'\\\\^'\\. 

In order to bound the last expression, define a Hamiltonian H(s) = H(s) — 7 (s)J. Then, it is clear that 
\l/(s) is an eigenvector of H(s) with eigenvalue and all other eigenvalues are at least A in absolute value. 
Therefore, according to Lemma 1^21 

m ^\\HJ^ \\H'\\+ max, ^ < 2\\H'\\ 



X ~ A ~ A 

and we obtain 

7 "(0) < \\H"\\+A\\H'\\ 2 /X 
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